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We review a recent approach for the simulation of many-body interacting systems based on 
an efficient generalization of the Lanczos method for Quantum Monte Carlo simulations. This 
technique allows to perform systematic coirections to a given variational wavefunction, that 
allow to estimate exact energies and correlation functions, whenever the starting variational 
wavefunction is a qualitatively correct description of the ground state. The stability of the 
variational wavefunction against possible phases , not described at the variational level can be 
tested by using the "effective Hamiltonian" approach. In fact Monte Carlo methods, such as 
the "fixed node approximation" and the present "generalized Lanczos technique" (Phys. Rev. 
B 64,024512, 2001) allow to obtain exact ground state properties of an effective Hamiltonian, 
chosen to be as close as possible to the exact Hamiltonian. thus yielding the most reasonable 
estimates of correlation functions. We also describe a simplified one-parameter scheme that 
improve substantially the efficiency of the generalized Lanczos method. This is tested on the 
t — J model, with a special effort to obtain accurate pairing correlations, and provide a possible 
non-phonon mechanism for High temperature superconductivity. 

1 Introduction 

Despite the tremendous progress of computer performances the general task of determining 
the ground state wavefunction of a many-electron system is still far from being settled For 
instance, even for simplified models on a lattice, there is no general consensus on the 
ground state properties of a system of about 100 electrons on L 100 sites. The most 
striking example is the so called t — J model: This model is still a subject of intense 
numerical studies, due to its possible relevance for High Tc superconductivity^^^. The 
Hamiltonian reads: 



H = J 2_^lSi-Sj - -n,nj j - t 2^ c[^Cj>, (1) 

where — cl ^ {1 — fii^g), hi = X^o- "■'■o- '^^e electron density on site i. Si = 

Scrcr' cr'''o->cr'^j>o-' ^V^^ Operator and r^.cr' are Pauli matrices. In the following 

we consider N electrons on L sites, with periodic boundary conditions, (PBC), in order to 
minimize size effects. 

After many years of intense numerical and theoretical efforts there is no general con- 
sensus on the properties of this simple Hamiltonian and of the related Hubbard model. In 
particular according to density matrix renormalization group (DMRG) studies d-wave 
superconductivity is not stable in this model, whereas a ground state non uniform in den- 
sity (with so called "stripes") is found. Several QMC studies provide controversial results, 
most of them indicating a superconducting behavior, and some of them^, indicating the 
opposite. 
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The reason of the above controversy, can be easily explained within the straightforward 
variational approach. Whenever a model Hamiltonian cannot be solved exactly either nu- 
merically (with no sign problem) or analytically, the most general and reasonable approach 
is an approximate minimization of the energy within a particular class of wavefunctions, 
for instance also DMRG can be considered a variational approach with a particularly com- 
plicated variational wavefunction obtained by DMRG iterations. However, within the vari- 
ational approach, one faces the following problem: for large system size L the gap to the 
first excited state scales generally to zero quite rapidly with L. Thus between the ground 
state energy and the variational energy there maybe a very large number of states with 
completely different correlation functions. In this way one can generally obtain different 
variational wavefunctions with almost similar energy per site, but with completely differ- 
ent correlation functions. It is easily understood that, within a straightforward variational 
technique, there is no hope to obtain sensible results for large system size, unless for for 
system with a finite gap to all excitations, such as spin Uquid^, or band insulators.. 

In the following we are trying to argue that a possible solution to the previous limitation 
of the variational technique is provided by what we call in the following "the effective 
Hamiltonian approach". 

This approach relies on the following assumption: 

"Among similar Hamiltonians with local interactions the ground state correlation func- 
tions depend weakly on the details of the Hamiltonian, in a sense that similar Hamiltonians 
should provide similar correlation functions"'. In this way the ground state of an effective 
Hamiltonian (such as the fixed node Hamiltonian*) that can be solved exactly by Quantum 
Monte Carlo schemes can be used as a variational state of the desired Hamiltonian, in this 
way providing not only a good variational energy but the most reasonable estimate of cor- 
relation functions, as long as the variational energy obtained is close -but not terribly close 
as in the straightforward variational approach- to the exact ground state energy. 

The paper is based therefore on the recent numerical advances for solving approxi- 
mately model Hamiltonians on a lattice: the fixed node*, and the "generalized Lanczos 
technique"^, that allows to improve systematically the variational energy provided by the 
effective Hamiltonian approach, by combining in an efficient way the power of the Lanc- 
zos variational technique with the "effective Hamiltonian approach". Trough all the paper 
and pictures we will use "FN" to indicate the "fixed node approach", , whereas "SR "will 
indicate the . "stochastic reconfiguration method" used to apply the "generalized Lanczos" 
scheme. In the first part we describe the Lanczos technique, then we derive the effec- 
tive Hamiltonian approach in a slightly more general way than the standard "fixed node" 
method. Finally we show that the mentioned "generalized Lanczos method" represents a 
very efficient implementation of both the previous techniques- Lanczos and fixed node- on 
a lattice. We also point out some slight but important improvements and simplifications 
to the most recent formulation of the "generalized Lanczos scheme"^. In the last section 
before the conclusion we show some example on the t-J model, where the "effective Hamil- 
tonian approach" is clearly useful, as the pairing correlation functions appear to be rather 
independent from the initial variational guess, even for large system size L ~ 50 and small 
J/t. 
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2 The Lanczos technique 



The Lanczos technique represents a remarkable improvement of the power method used to 
filter out systematically the ground state component of a given initial wavefunction by 
an iterative technique. The power method is based on the following equation: 

Vo> ^ (A/ - /^nVc) (2) 

where A is a suitable large shift to ensure convergence to the ground state for large p, I 
is the identity matrix and |V'o) the ground state of if. At a given iteration p, after apply- 
ing just p powers of the Hamiltonian, a much better wavefunction tjjp can be obtained by 
combining, with proper coefficients a^, the states obtained with the power method in the 
previous iterations: 

p 

iV'p) = (i+E"'^'f^')i^G> (3) 

fe=l 

with parameters {ak} for k — I,-- - ,p minimizing the energy expectation value 
{iljp\H\il)p) / {ipp\tpp) . For any p it is simple to show that the wavefunction (3) corre- 
sponds exactly to apply p Lanczos step iterations to the initial wavefunction [iha). The 
iJ— polynomial of degree p which is applied to the initial state ^po, can be generally fac- 
torized in terms of its roots Zi'. 

[l + j2'^kH'')=flil-H/z,) (4) 
fe=i 1=1 

This decomposition will be particular important for applying statistically the Lanczos tech- 
nique with the Stochastic Reconfiguration (see later). As it is clear from Fig. (1), the Lanc- 
zos method converges very quickly to the ground state wavefunction especially when a 
particularly good "guess" is used for iJjq- 

Whenever the ground state wavefunction is approached | {tpo | V'p ) P / ( V'p I ) ^ = 1 — ep> 
with Cp — > for larger p, with the energy approaching the exact value with corrections 
~ Cp. On the other hand, the variance ap of the Hamiltonian on the approximate state ipp 

al = {i^AH^\xPp) - {^Ppmpf = 0{ep) 

is going to zero in the limit when ipp is the exact eigenstate tpo with the same corrections 
proportional to e^.: It is clear therefore that a very stable evaluation of the energy can be 
done by using few Lanczos steps values of the energy and the corresponding variance. 
Then, by performing simple extrapolation (hnear or even polynomial), the exact ground 
state result is easily estimated provided the energy-variance values are close to the Un- 
ear regime (see Fig.l). The same scheme can be applied even for correlation functions^, 
and represents one of the most simple and effective methods to estimate exact correla- 
tion functions with few Lanczos steps (i.e. with a minor computational effort) whenever 
the variational wavefunction tpc is particularly good, i.e. is close to the linear energy vs. 
variance regime. Such property of the variational wavefunction can be satisfied even for 
system size L ~ 100^. 

The initial wavefunction to which the Lanczos and the following techniques will be 
appUed can be written as follows ^°: 

\iPG) = \^p=o) = PoPNJ\D). (5) 
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Figure 1. Energy < > vs. variance < > — < H >^ of the Lanczos technique for different initial 
wavefunction tpG- Here n represents the number of iterations. Lower variance is always obtained for larger n. 
The zero variance limit is the exact results. 



where \D) is a BCS wavefunction, which is an exact eigenstate of the following Hamilto- 
nian: 

Hbcs = Ho + (At + A) (6) 

where Ha — J2 c|. „c.k,cr is the free electron tight binding nearest-neighbor Hamilto- 

k.{7 

nian, = ~2{coskj. + cosky) — /i, /i is the free-electron chemical potential and At 
creates all possible singlet bonds with d-wave symmetry being Mij, Mij not restricted 
to nearest neighbors, but exhaustively parametrized with a reasonable number of varia- 
tional parameters as described in ^. and Pq are the projectors over the subspaces with 
a fixed number N of particles and no doubly occupied states. Finally the Jastrow factor 
J — exp (1/2 J2i j ~ ) couples the holes via the density operators fii and con- 
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tains other few variational parameters. We note here that by performing a particle-hole 
transformation on the spin down ct — > (— l)'Cj,|, the ground state of the BCS Hamil- 
tonian is just a Slater-determinant with N = L particles This is the reason why this 
variational wavefunction can be considered of the generic Jastrow-Slater form, a standard 
variational wavefunction used in QMC. All the mentioned variational parameters are ob- 
tained by minimizing the energy expectation value of H over "ipG^- 

Using the particle-hole transformation, it is also possible to control exactly the spurious 
finite system divergences related to the nodes of the d-wave order parameter. 



3 The effective Hamiltonian approach 



In a discrete Hilbert space defined for instance by configurations x of electrons with defi- 
nite positions and spins we consider any Hamiltonian H with real matrix elements H^'^x 
and any real wavefunction iI)g{x) assumed to be non zero for each configuration x. 

By means of the wavefunction i/^g -hereafter called the guiding wavefunction- we can 
define a two parameter class of Hamiltonians -f?^^ depending on 7 and r: 

( Hx,x + {l + l)Vsi{x)+r{l+-f)eL{x) for x' = x 

Hlj^=l Hx',x iix' ^x and il)G{x')Hx',x/i^G{x) <0 

\ - 7-ffx',x if x' ^ X and 'iI)g{x')Hx',x/'4>g{x) > 

(8) 

where the local energy ei(x) is defined by: 

^l{x) = ^ iJG{x')Hx',x/^G{x) (9) 

x' 

and the so called sign-flip term Vgf {x) introduced in * is given by considering the sum of 
all the positive off-diagonal matrix elements appearing in the local energy. The effective 
Hamiltonian has the same matrix elements of the Hamiltonian H for all off-diagonal 
matrix elements that do not frustrate the guiding function signs, the other ones are taken 
into account by proper modification of the diagonal term. 

The following properties are almost an immediate consequence of the above defini- 
tions: 



i) for7=-l H = H}^,. 

ii) for r = — 1/(1 -|- 7) and 7 7^ —1 the ground state of -f^^jv '•^^ guiding wavefunction 
itself with zero ground state energy, namely Hp^\tpG) = 0. 

iii) if = if^^-(l + 7)^ 

iv) El{x) = Y.X' ^Jg{x')HIn/iJjg{x)= eL{x){l + r(l + 7)) where ELix) is the local 
energy of the effective Hamiltonian Hpj^, whereas eL{x) = J2x' '4^g{x')H/'iPg{x)- the 
corresponding one for H. Moreover: 

(v) for 7 > the ground state ijj^^ [x) of i?^^ may be chosen to have the same signs 
of the guiding wavefunction, namely tpG{x)^FN{x) > for any configuration x. This 
follows by doing a unitary transformation of the basis \x >= Sign [^G(x)]|a; >, in which 
the off-diagonal matrix elements of the Hamiltonian Hpjy ^/ ^ < are non-positive. Thus 



5 



the Perron-Frobenius theorem holds implying that a ground state wavefunction (in prin- 
ciple there maybe degeneracy) can be chosen to satisfy ipQ^(x) > in the new basis , 
which finally proves (v) in the original basis. The statement (v) suggests that the effective 
Hamiltonian H]^^ represents the lattice counterpart of the fixed node (FN) hamiltonian, a 
well known approximation for continuous models.^ Furthermore, provided the matrix ele- 
ments of the hamiltonian H or H^^ satisfy an ergodicity property (namely that any two 
arbitrary configurations x and x' can be always connected by a suitable large number M of 
hamiltonian powers {x'\H'^\x) ^ 0), then a more restrictive property holds: the ground 
state is unique for any 7 > 0. This implies immediately that: 

(vi) the ground state energy E{"f) of the fixed node hamiltonian -ff^jv ™ analytic func- 
tion of 7, due to the finite size gap separating the unique ground state from the first excited 
state. We assume in the following that this very general property holds for the given hamil- 
tonian a condition which is not restrictive, also considering that if ergodicity is not satisfied, 
all previous and the following considerations hold in all the subspaces of configurations x 
ergodically connected by the powers of the hamiltonian. 



By using Green Function Monte Carlo the ground state energy £'(7) can be very 
efficiently computed for 7 > as all the matrix elements of the importance sampled 
Green function G^,^ = %Pg{x') [Mx',x - {Hfn)x',x\ /^g{x) are all positive for large 
enough constant shift A. This is obtained by averaging the local energy < El{x) > over 
the configurations x generated statistically by the Green function G^^ with a standard 
algorithm.^'^^'^'^ Notice also that, by property (iv), the local energy El of this fixed node 
hamiltonian is proportional to the local energy of H and therefore this computation 
satisfy the so called zero variance property: both E^ and ei, have zero statistical variance 
if ipG is an exact eigenstate of H. 

For r = Hj,^ reduces to the standard fixed node hamiltonian defined in ^ ( 7 = 0) 
and extended to 7 7^ in Thus a rigorous theorem holds relating the ground state 
energy £^(7) of the fixed node ground state i^pj^ of Hp^^, to its variational expectation 
value £^^(7) = {tppj^\H\ippj^) on the hamiltonian il: 

£^^(7) < ^^(7) < {MH\^g) (10) 

Using property (i) we therefore notice that by increasing the value of r from the vari- 
ational value r = —1/(1-1-7) up to r = the ground state of the fixed node hamiltonian 
Hpj^ becomes a variational state with lower energy expectation value. This implies imme- 
diately that the fixed node effective hamiltonian is more appropriate to describe the ground 
state of H. 

In the continuous case r cannot be extended to positive values because the local energy 
bl may assume arbitrary large negative values close to the nodes, and the best variational 
energy can be actually obtained just for r = (since for r = the fixed node gives the 
lowest possible energy compatible with the nodes of the guiding function). In a lattice case 
such a theorem is missing, and there is no reason to expect that r = is just the optimal 
value. 

A simple and efficient scheme to compute a variational upper bound of the energy for 
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any r is described in the following paragraphs. Using property (iii) 

EMl) = {^InW-^ - (1 + 7)^IV'J^) = E{^) - (1+7)^ (11) 

where in the latter equality the Hellmann-Feynmann theorem has been used. By using that 
il^jY depends linearly on 7, the well known convexity property of E{'y) holds^^ : 

Therefore the expectation value Efn{i) of the hamiltonian H on the fixed node state is a 
monotonically increasing function of 7), as clearly = — (1 + > 0. The 

best variational estimate is obtained therefore for 7 = 0, as in the conventional scheme. 

The extension to finite 7 is however convenient to provide better variational estimates 
of Ej,~^, which in fact maybe sizable lower than the standard estimate Epr^{0) < E{0) 
for r = 0. This extension allows also to make a rigorous upper bound of Epj^ also in the 
case r > 0, without missing the zero variance property. In fact, always by the convexity 
property of £^(7), 

dEij) E{j) - E{0) 

"^'^=°- J (1^) 

we finally get that at the best variational condition 7 = 

Efn{0) < E{0) - (£(7) - E{0))/-f. (14) 

For r = the above upper bound improves also the previously known value (10), at least 
for 7 small enough where the above inequality becomes a strict equality. 

In practice, since the energy as a function of 7 is almost linear a very good estimate 
can be obtained using the above inequality even for 7 = 1, as shown in Fig. (2) for a 
test example on the i — J model, where it is also clear that the variational energy can be 
improved by turning on the parameter r. 



4 The generalized Lanczos 

The optimization of the parameter r is rather problematic within the scheme of the previ- 
ous section especially when few Lanczos steps are applied to the guiding function and the 
dependence of the energy as a function of r cannot be resolved within available statistical 
errors. Though the energy maybe rather insensitive to r, the behavior of correlation func- 
tions, may strongly depend on it, especially when the guiding function shows some insta- 
bility towards different phases not described at the variational level. Within this approach 
the instability of the guiding function is characterized by the existence of a considerable 
number of configurations x with local energy (x) much below the average and with cor- 
relation properties much different than the average. By increasing r these configurations 
will have larger and larger weight in the fixed node ground state V'fjv (since they have 
much lower-energy diagonal term) and wiU display clearly the possible instabilities of the 
variational wavefunction tpc- 

The sign-flip term Vsf{x) is divergent whenever the guiding function is exceedingly 
smaU (i.e. close to the nodes or finite-size lattice pseudo-nodes of ipa), thus requiring 
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Figure 2. Variational energy of the t-J hamiltonian as a function of the pai'ameters r and 7, for the BCS-guiding 
function (5), without any Lanczos improvement. The 7 — > limit in the right panel con'esponds to the expectation 
value {■y) = ("i/'^jy I^^IV'J'jv) for 7 = where V'Jjv ground state of the effective hamiltonian Hp^^. 
Each point, due to inequality (14), represents an upper bound for {-y = 0) and, clearly, for the ground state 
of H. All the estimates reported here are much better than the standard r = lattice fixed node upper bound 
E{f = 0)8 for £^^(7 = 0): £(7 = 0) = -0.77580(2) much above the upper energy scale. The value (SR) 
obtained with the "generalized Lanczos" described in the following sections is also shown for comparison. 



an infinite shift A^"', because for the statistical implementation of the power method the 
diagonal term A - {Hj^j^)^,^ = A - H^.^ - (1 + 'y)Vsf{x) - r(l + j)eL{x) (see Eq.8) 
has to be non negative. For r = — 1/(1 + 7), in the variational case, a better approach, 
but similar in spirit, is obtained by sampling^^ the square of the variational wavefunction 
ipQ with a different Green function. This following importance sampled Green function is 
used for the statistical implementation of the power method: 



^^'■^ = I ~^i^GW){Hl^)l.. J^G{x) for x'^x 
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where Zx is a normalization factor obtained by setting z^'Gl, ^ = Zx, namely:. 

Zx = K-eL{x) + {\ + 7) V,/ {x) (16) 
In this way it is straightforward to show that: 

Y,Gl,J^l^G{x)\^ = \'4^G{x)\^ (17) 

X 

Thus the importance- sampled Green function C' maybe used to generate configurations 
that sample the variational wavefunction square. The advantage of the present approach is 
evident since the diagonal term of the Green function does not contain the sign-flip term, 
and a finite reasonable A can be used. For instance in the t — J model A can be set to 
zero. Instead a zero shift is not allowed for the importance sampled Green function of the 
effective hamiltonian itself: 

Gfn = ^g{x') [A - {Hl^)x>.x] /^g{x) (18) 

which performs the same task for r = —1/(1 + 7), but with a less efficient infinite A 
scheme^^. 

In the following, within the spirit of the "effective hamiltonian approach", the varia- 
tional wavefunction is improved by tuning a parameter r proportional to the local energy, 
in order to modify and improve the effective hamiltonian H^.^^, whose ground state is just 
V'G for f = ^1/(1 + 7)- This parameter is then changed in order to be as close as possible 
to the true hamiltonian for 7 > 0, when computations free of sign problem are possible. 
Indeed in order to improve -ff^jv it is very useful to notice that H'^jf = H, the exact hamil- 
tonian, for 7 = — 1 and any non-zero r. Thus at finite positive 7 an optimal variational 
parameter r can be used, that on a lattice, maybe significantly different from the fixed node 
value r = 0, since this value represents the optimal one only in a continuous model, when 
there exists a rigorous proof that r = provides the minimum possible energy. 

In order to determine a feasible scheme for the optimization of r in the lattice case, we 
need to implement small modifications of the Green function (15). We notice that there are 
two important changes of this Green function that are easily implemented: 

4.1 One lanczos step improvement 

In this case the Green function (15) is modified by: 

GJls = r.'GlJrx (19) 

where Vx = 1 + acL (x) . After applying statistically the above Green function, after a large 
number of iterations the configurations x, will be distributed according to the weight (not 
necessarily positive): 

tjjG{x)tl)i {x) 

where 

= (1 + aH)\^G) =Y^rxipG{x)\x) (20) 

X 

is the first Lanczos step wavefunction as described in Eq. (1). Since the Lanczos iteration 
improves the wavefunction and the factor has not a definite sign on each configuration 
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X, it is clear that the phases of the ground state wavefunction are much better represented 
by the signs of rxtpcix) rather than by the ones corresponding to tpcix). The parameter 
a = ai/ao can be determined by satisfying the SR conditions^: 

{iPG\H{ao + aieL)|Vn) = {i^GH{A - 

(VgIK + aiBLMn) = (V-gKA - H)\i,n) (21) 

where ai,i = 0, 1 are computed statistically at any given iteration n in order to improve 
the SR state Tx'tpnix), until convergence is reached for large n. In this case tpn{x) is 
independent of n and statistically equal to ipG, whereas a will converge (statistically) to the 
exact one Lanczos step value. Once this value is determined the energy expectation value 
over tpi can be evaluated by statistically averaging the local energy eL{x) corresponding 
to tpG (and not to ipi), providing a substantial reduction of computational effort. In this 
case, since the value of 7 is immaterial for the statistical averages, it is more convenient to 
use 7 = 1, that minimizes statistical fluctuations. 

In general, the use of the SR conditions^ allows to obtain the energy and correlation 
expectation values of the Lanczos step wavefunction i/jp, by using a guiding function 
ipG containing only p — 1 powers of the Hamiltonian, e.g. |V'g) iV'p-i)- The use 
of IV'p-i) as a guiding function for sampling rpp may not be the optimal choice. In the 
following we describe a guiding function with better nodes than ipp-i but with the same 
number p — 1 of hamiltonian powers, that will be used in the following sections whenever 
the method SR will be applied. 

Using the root decomposition (4) of the H — polynomial defining the p— Lanczos step 
wavefunction |^p), we can single out any real root Zk and similarly to the first Lanczos 
step case: 

il)p{x) = r^'^Gix) with 

r-x = 1 - eL{x)/zk 
\^G)^l[{l-H/zi)\^G) (22) 

The new local energy eL{x), obtained with the new guiding function, will keep into ac- 
count the phases of the p— Lanczos step wavefunction exactly. In this way, within this 
decomposition, it is clear that the best guiding function i/jg of the previous form, is ob- 
tained by choosing the real root Zk such that: 

<l-eL(x)/zfe> (23) 

is as far as possible (on average over i/jg) from the zero value. This condition (23) will min- 
imize the sign changes of ipcix) to obtain tpp{x) = (1 — e^ix)/ Zk)i^G{x), thus providing 
the best possible phases that we can safely obtain with p — 1 powers of the hamiltonian 
appUed to the bare V'g- 

4,2 Fixed node improvement 

In this case the Green function is modified similarly: 

G'fn = r,.Gl,JSgn{r,) (24) 
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It is easily obtained that for r^; = 1 — g^^'^.^ ^j^^j large shift A, the effective 

hamiltonian -ff^^ (8) is indeed considered, as for A ^ 00 the matrix elements of Gfn 
(18) coincide with the ones defined above for AG'pj^. up to 0( j). 

In particular for r = 0, and 7 = we recover the standard fixed node*. Notice also 
that, if the hamiltonian is free of sign problem Vsf{x) = and the fixed node is exact. 
Then the choice r = provides the exact sampling of the ground state of H even for finite 
A, as the factor is proportional to Zx (16) and simplifies in (18,15). 

4.3 Generalized Lanczos 

Using the above Green function (24), the parameter r = "^^"j^^""^"^ , a single param- 
eter at any order p of the Lanczos iterations, is optimized using the SR conditions(21) 
with il^n now depending explicitly on n and differing from the initial guiding function 

tpG- Txi^n{x) = {G'p^YipG- These conditions provide, as mentioned before, ao, 
statistically.^: However, in this case, the parameter r, determined by the SR condition, 
may not coincide with the lowest possible energy condition. A further modification of the 
Green function'^ 

G'r, = rx.Gl,J\rx\'~^'^Sgn{rx) (25) 

that interpolates between the Lanczos limit (19) for = (when the SR conditions coin- 
cide with the Euler condition of minimum energy) and the fixed node limit (24) for t] = 1 
allows to overcome this difficulty, as we get closer but not exactly equal to the Lanczos 
limit, and one can obtain even lower variational energies.^ 

For the t — J model we avoid to consider here this extra-complication, since the SR 
conditions (21) have been tested to coincide almost exactly with the Euler conditions of 
minimum energy (see Fig. 2) even for 77 = 1 at least for A = 0. As shown in the same figure 
the SR may also provide a slightly lower energy than the corresponding one obtained by the 
best r effective hamiltonian Hppf, because for small A the factor Vx in Eq.(24) may change 
sign and can correct also the phases of the wavefunction and not only the amplitudes. This 
is also the reason to work with the minimum possible shift A. In principle it is possible 
to further improve the variational energy and the nodes of the sampled wavefunction, by 
performing the reconfiguration scheme each kp steps, with an effective Green function: 

G'k,=rx'iG'')''jJSgn{rx) (26) 

For 7 = 1, it is possible to work with kp > 1 and with reasonable statistical fluctua- 
tions (that increase obviously with kp). By increasing kp the factor provides non triv- 
ial changes to the phase of the wavefunction with corresponding improvement in energy 
expectation value. We have not systematically studied this possible modification of the 
method so far. This extension to fcp > 1 should be clearly useful for model hamiltonians, 
such as the Hubbard model at strong coupling, when a large shift A is required for the 
convergence of the method. 

For A = or finite, the coefficient r in the factor may have little to do with the 
coefficient appearing in Hj,^, but, even at finite A, an effective hamiltonian can be still 
defined"^, which is qualitatively similar to Hp^. In the following discussions we will not 
consider the difference between the finite A effective hamiltonian and the infinite A one 
(8) because it is irrelevant for our purposes. 
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At each iteration p of the generahzed Lanczos the special guiding function described 
in Eq. (22) is used, yielding optimal phases as close as possible to the p— Lanczos step 
wavefunction. As far as the remaining parameter 7, this is restricted to be positive for 
statistical reasons (no sign problem). Clearly from property (12), the smaller is 7, the better 
is the variational energy but increased fluctuations occurs for computing the SR conditions 
(21). On the other hand, the Green-function shift A has to be taken as small as possible, 
compatibly with A — H^^x > for any x. in order to further improve the efficiency of 
the power method. Within the SR method by minimizing at best the parameters 7 and 
A (or increasing kp) we can further improve this technique, in a practical scheme. The 
optimization of the parameter r, since it affects a change in the effective hamiltonian is 
particularly important for correlation functions. Instead all the other parameters (including 
rj or kp for instance) may help to obtain slightly lower variational energies, but are in 
general much less important. The variational SR results for the i — J model, described in 
the following sections, are obtained with 7=1/4 and A = and refer to the fixed node 
Green function (24), whereas the symbol FN will always refer to the standard fixed-node 
case A ^ cxD, 7 = r = 0. 

5 Results on the t- J model 

We consider the pairing correlations in the f — J model for square clusters with periodic 
boundary conditions: 

= 4,Al + (27) 

a] j creates a singlet pair in the sites i, j. On each lattice we take the first nearest neigh- 
bor pair i,j fixed and move k, I parallel or perpendicular to the direction In all cases 
studied the parallel correlations are positive and the perpendicular ones are negative, con- 
sistent with a (i— wave symmetry of the pairing. The existence of phase coherence in the 
thermodynamic limit is obtained whenever Pij;k,i remains finite for large distance sepa- 
ration between the pair i, j and k, I. A systematic study has been reported in^^. Here we 
focus only on few test cases to show the power of the method, and the importance to work 
with an effective hamiltonian Hpj^ with a single variational parameter r as described in 
the previous section. For all cluster used the distance between pair i, j and pair k, I refers 
to the minimum one between \Ri — Rk\, \Ri — -R; |' \Pj ^ Rk \ and \Rj — Ri\. Only for the 
6x6 we use the so called Manhattan distance = |a;| + since the pair {k, I) in this 

case is moved in both perpendicular directions. First the pair (fc, /) is translated paraUel to 
the a:;-axis up to the maximum distance allowed by PBC, and then (for the 6x6) the pair 
(/c, I) is moved parallel to the y— axis. 

First of all, whenever the initial variational wavefunction used is qualitatively correct 
(5), few Lanczos iterations are really enough to obtain exact ground state properties. This 
is clearly shown in Fig. (3) where the exact results coincide within few error bars with the 
variance extrapolated results, that in turn are very close to the p = 2 Lanczos wavefunction 
results. However for larger system when the solution is not known, few Lanczos iterations, 
though systematically improving the energy, cannot change qualitatively the pairing corre- 
lations of the initial wavefunction, and in general the variational approach is not reUable. 
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Figure 3. Pairing correlations in tiie 26 lattice for 4 holes in the J/t = 0.5 t~ J model for the variational Lanczos 
technique as compared with the exact result obtained with exact diagonalization. The variance extrapolated values 
are obtained using only the p = 0, 1, 2 results available with the statistical algorithm also for much larger system 
size. 



In order to show this effect, we have used two different variational wavefunctions on 
a 6 X 6 4-holes J/t = 0.5 cluster, and improved both initializations with the methods 
described in the previous section: the pure variational Lanczos technique, the standard 
fixed node (FN) and the "generalized Lanczos method" (SR), within the simplified scheme 
considered before. For one wavefunction initialization, the BCS variational parameters 
are optimized by minimizing the energy, for the other one we have reduced to a very small 
value ~ 10^'' the corresponding variational parameter A bcs in (6), just in order to remove 
the degeneracy of the free-electron determinant in the 6x6. This choice yields a variational 
wavefunction with definite quantum numbers and with small pairing correlations. 

We see in Fig. (4), top panels, that the Lanczos technique is very much dependent on 
the two different initial choices, even though the energy is in both cases very much im- 
proved by few Lanczos iterations. As shown in Fig. (5), the variance extrapolated results 
of the energy are consistent for both initial wavefunctions. On the other hand the pairing 
correlations remain inconsistent for about a factor two at large distance. 

In this example we clearly see the limitation of the straightforward variational tech- 
nique: within a very similar energy (e.g. the extrapolated ones) the pairing correlations 
maybe even qualitatively different. 

A completely different behavior is obtained as soon as the FN is applied (middle panels 
in Fig. 4). The energy improvement within this technique is apparently marginal compared 
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Figure 4. Pairing correlations in the 6x6 lattice for 4 lioles in the J/t = 0.5 t — J model. Left panels and 
right panels refer to two different initial guiding functions with or with vanishing small d-wave order parameter 
respectively. The latter is used in order to remove the degeneracy of the free electron Slater-determinant. The 
panels at different raws refer to different methods, as a function p of the hamiltonian powers used to evaluate 
the local energy ej^, required by all the methods: the larger is p. the more ( for p > 2) computationally 
demanding is the calculation. The VMC values (red trangles) are plotted in all panels for comparison. 



to the standard Lanczos technique (see Fig. 5). Instead the behavior of pairing correlations 
is much better, and akeady the simple fixed node approximation applied to the pairing 
correlations is rather independent of the initial wavefunction. The only drawback of this 
technique is that when systematic improvements to the variational wavefunction are imple- 
mented ( larger p in the figure), the convergence properties are not behaving so accurately, 
as one could expect from the convergence of the energy reported in Fig. (5). In partic- 
ular, even at the most accurate level of this fixed-node approximation -namely the fixed 
node over the two Lanczos step wavefunction- the two different initializations give pair- 
ing correlations differing by about 20% at the largest distance. This is much better than 
the straightforward Lanczos variational technique (this difference was about 70% for the 
corresponding two Lanczos step wavefunctions) but is not satisfactory enough. 

The reason of such behavior is easily understood in terms of the effective hamiltonian 
approach. In a lattice case it appears really important for correlation functions to optimize 
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Figure 5. Variational energies obtained witli various metliods as a function of tlie variance per site jl? of the 
p— Lanczos step wavefunction (VMC), wliicli is improved either with standard fixed node (FN) or the generalized 
Lanczos (SR), with the simplified and efficient scheme described in the previous section. The values at zero 
variance are exQ'apolations with a quadratic least square fit. 



the parameter r appearing in the effective hamiltonian (8) and not just taking the FN ansatz 
r = 0. This optimization scheme is particularly important whenever some correlations 
that are not included at the variational level (or much depressed as in the case studied) 
are increasing as we go down in energy with the help of the improved p — 1 {p > 1) 
Lanczos step guiding function. In general for larger p the parameter r increases, thus the 
SR scheme provides correlation functions substantially different and more accurate than 
the FN. In the bottom panels it is remarkable that, after applying only 3— steps of the 
SR technique, both initializations provide the same results within error bars (< 3%j at 
the largest distance. These results can be considered benchmark accurate calculations of 
pairing correlations in the QxQ cluster These pairing correlations clearly indicate a robust 
wave superconducting ground state in the t — J model, at least for this J /t ratio. In this 
example we notice that correlation functions, in the effective hamiltonian approach, begin 
to be consistent within 5% whenever the variational energy is accurate within ~ 1%, that 
is at least one order of magnitude better than a straightforward variational technique like 
the Lanczos one. 

Of course for larger size, consistent correlation functions, i.e. independent from the 
initial wavefunction with or without Abcs^ can be obtained for a larger number p of SR- 
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Figure 6. Pairing correlations in tlie 50 site lattice for 8 holes in the J/t = O.l t — J model. Left panels and 
right panels refers to different initial guiding function with or without d-wave order parameter respectively. The 
pairing coiTelations for both calculations are consistently small at the most accurate level of approximation (SR 
p = 3). 

iterations. Here we report a sample case for a 50 site cluster at small J/t = 0.1. We see 
in Fig. (6) that the sizable pairing correlations present in the variational wavefunction with 
Abcs 7^ 0, represents just an artifact of the variational calculation. At the third step, of 
the SR technique, when, as shown in Fig. (7) we reach an accuracy in energy below 1% 
(assuming that the variance extrapolated energies-both consistent- are exact), the pairing 
correlations are again consistent within few error bars, and clearly vanishingly small at 
large distance. 

6 Conclusions 

We have shown that within a brute force variational technique, such as the Lanczos method 
for few iterations, it is hard to obtain accurate values of correlation functions unless the en- 
ergy accuracy is far from the present possibilities, at least in two dimensions. An accuracy 
of about one part over 10* in the energy would be probably possible with at least 10 
Lanczos steps or 100000 states in DMRG 2D calculations for systems of about 100 sites 
with periodic boundary conditions. This kind of accuracy maybe enough to obtain consis- 
tent correlation functions even within these two variational methods, but is far from being 
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Figure 7. Variational energy as a function of the variance per site jl? for tlie p— Lanczos step wavefunction 
(VMC), which is improved by the "generalized Lanczos method" (SR). The best variational SR p = 3 energies 
are indicated by the arrows. 



possible at present. 

We have shown instead that a qualitatively different and very promising approach, 
based on the optimization of an effective hamiltonian, rather than adding more variational 
parameters in a brute force variational scheme, appears to be very useful to control corre- 
lation functions. The idea is based on the "effective hamiltonian approach" described in 
the introduction. In this scheme it is assumed that between similar Hamiltonians, the corre- 
lation functions of their ground states should be also similar. The SR technique, allows to 
systematically improve the effective hamiltonian considered even compared to the lattice 
fixed node one®, with an iterative scheme very similar to the Lanczos one, thus the name 
"generalized Lanczos". 

Within this scheme it is clear that there are robust pairing correlations in the t — J 
model at sizable but not unphysical value of J/t^'^ . However there exists a critical value 
{J/t)c > 0.1 below which pairing correlations are clearly suppressed. The existence 
of such a critical {J/t)c is clearly understood because at J/t = 0, the ferromagnetic 
instability takes place even at large doping 



17 



Acknowledgments 

We are indebted to F. Becca, L. Capriotti, A. Parola and E. Dagotto for many useful dis- 
cussions. This work was partially supported by MIUR, (COFIN-2001) and INFM (Pais- 
Malodi). 

References 

1. G. Baskaran, Z. Zou, and P.W. Anderson, Solid State Comm. 63, 973-975 (1987). 

2. F.C. Zhang and T.M. Rice, Phys. Rev. B 37, 3759,3763 (1988). 

3. S. Sorella, Generalized Lanczos algorithm for variational quantum Monte Carlo, 
Phys.Rev. B 64, 024512 1-16 (2001). 

4. S.R. White and D. Scalapino, Phys. Rev. Lett. 80, 1272-1275 (1998); ibidem Phys. 
Rev. B60, 753 (1999). 

5. C.Y. Shih, Y.C. Chen, and H. Q. Lin Phys. Rev. Lett. 81, 1294-1297 (1998). 

6. L. Capriotti, F. Becca, A. Parola and S. Sorella Resonating Valence Bond Wavefunc- 
tions for Strongly Frustrated Spin Systems Phys. Rev. Lett. 87, 098201 (2001). 

7. N. Trivedi and D. M. Ceperley Phys. Rev. B 41, 4552 (1990). 

8. D. F. B. ten Haaf, J. M.J. van Leeuwen, W. van Saarloos, and D. M. Ceperley, Phys. 
Rev. B 51, 13039 (1995) 

9. see e.g. D. M. Ceperley J. Stat. Phys. 63, 1237 (1991). 

10. C. Gros, Phys. Rev. B 38, 931 (1988). 

1 1. H. Yokoyama and H. Shiba, J. Phys. Soc. Jpn. 57, 2482 (1988). 

12. K. Runge, Phys. Rev. B 45, 12292 (1992); ibidem , 7229 (). 

13. M. Calandra and S. Sorella Numerical study of the two dimensional Heisenberg model 
using a Green function Monte Carlo technique with a fixed number of walkers, Phys. 
Rev. B 57, 11446-11456(1998) 

14. S. Sorella and L. Capriotti, Phys. Rev. B 61, 2599 (2000). 

15. This follows by writing 7 = jr/i + (1 — p)"i2, for any 71, 72 and < p < 1, thus 
finding a variational lower bound for E{"f) > pE{ji) + (1 — p)E{j2) because the 
ground state energy £^(7) of H'' is certainly bounded by the minimum possible energy 
that can be obtained by each of the two terms in the RHS of the following equation: 
H'* = The aboveinequality represents justtheconvexityproperty 
of £;(7). 

16. C.S. HeUberg and E. Manousakis, Phys. Rev. B 61, 1 1787 (2000). 

17. S. Sorella, G. Martins, F. Becca, C. Gazza, A. Parola and E. Dagotto, Superconduc- 
tivity in the two dimensional t-J model cond-mat/01 10460 , (2001) 

18. F. Becca and S. SoieM&Nagaoka Ferromagnetism in the Two-Dimensional Infinite-U 
Hubbard Model Phys. Rev. Lett. 86, 3396-3399 (2001). 



18 



